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ABSTRACT 

Massive dense cores (MDCs) are the high-mass equivalent of the so-called dense cores in nearby star-forming regions. With typical 
sizes of 0.1 pc, they could be either single high-mass protostars, or proto-clusters on the verge to form a cluster of low-mass stars. 
Here we present high angular resolution continuum observations obtained with the IRAM Plateau de Bure interferometer at 1.3 and 
3.5 mm towards the six most massive and young (IR-quiet) dense cores in the Cygnus X complex. Located at only 1.7 kpc, the 
Cygnus X region offers the opportunity of reaching small enough scales (of the order of 1700 AU at 1.3 mm) to separate individual 
collapsing objects, and thus to observe and constrain the result of the fragmentation process. The cores are sub-fragmented with a 
total of 23 fragments inside 5 cores. This is an indication that the monolithic collapse view may not apply to these cores. However, 
the most compact MDC, CygX-N63, stands out as a single object, suggesting it could actually be a single massive protostar with 
an envelope mass as large as ~ 60 M©. The fragments in the other cores have sizes and separations similar to low-mass pre- stellar 
and proto-stellar condensations in nearby protoclusters, and are most probably self-gravitating objects (M > Myir) and are either 
pre-stellar condensations or Class young stellar objects. A total of 8 out of these 22 protostellar objects are found to be probable 
precursors of OB stars in addition to CygX-N63, with their envelope masses ranging from 6 to 23 M© inside a FWHM of 3000 AU. 
The level of fragmentation is therefore globally higher than in the turbulence regulated, monolithic collapse scenario, but it is also 
not as high as expected in a pure gravo-turbulent scenario where the distribution of mass is dominated by low-mass protostars/stars. 
Here, the fractions of the total MDC masses in the high-mass fragments are reaching values as high as 28, 44, and 100 % in CygX- 
N12, CygX-N53, and CygX-N63, respectively, much higher than what an IMF-like mass distribution would predict. The increase 
of the fragmentation efficiency as a function of average density in the MDCs is proposed to be due to the increasing importance of 
self-gravity leading to gravitational collapse at the scale of the dense cores. At the same time, the observed MDCs tend to fragment 
into a few massive protostars within their central regions. We are therefore probably witnessing here the primordial mass segregation 
of clusters in formation. The physical origin of the fragmentation into a few high-mass objects is not clear yet, and will be further 
investigated in the future by investigating the kinematics of the protostars, and by measuring the importance of magnetic support, to 
check the two possible scenarios: early competitive accretion or magnetic supression of fragmentation during the global collapse. 

Key words. Stars: formation - Stars: luminosity function, mass function - ISM: clouds 



1. Introduction 

The detailed physical processes at the origin of high-mass stars 
(here defined as ionizing OB stars; i.e. earlier than B3 or > 
8 Mo) are still to be fully recognized and understood. Whether 
they form through a scaled-up version of the low-mass star for- 
mation process or require a specific process is still an open issue. 
The formation process might be particular for two main reasons. 
First, the radiation pressure of the ionizing protostar could stop 
accretion when it reaches a mass of ~ lOM©. With a moder- 
ate accretion rate 1Q~^ Me/yr), the star could at best reach 
~ 30 M© (see |Yorke & Sonnhalter|2002l and references therein). 
Second, OB stars are found to form in clusters and their for- 
mation could originate from collective effects. Competitive ac- 
cretion is one of the scenarios which require cluster forming 
conditions since the accretion efficiencies and rates depend on 
the gravitational potential of the proto-cluster ( [Bonnell & Bate] 
2006). The radiation pressure may not be a problem after all 
if the accretion rates are in practice much higher than usually 
measured in low-mass protostars (typically between 10"^ and 
10"^ Mo/yr). High infall rates could actually be expected in the 



classical theory of star formation applied to low-mass stars if the 
initial conditions for collapse significantly depart from a singu- 
lar r"^ sphere or if the collapsing conditions are reached in a 
dynamical way (see Larson 1969; Whitworth & Sum mers|1985[ 
Foster & Chevalier 1993 ; Henriksen et al. 1997 ; Chabri er|2003^ 
Alternatively, Mc Kee & Tan, ( ,20 02) propose that even in a slow, 
quasi-static evolution of star-forming clumps, large accretion 
rates could be due to the high level of turbulence. However, a 
continuous injection of turbulence in these high density regions 
would then be required. Turbulence is also usually believed to 
be at the origin of the hierarchical structure of molecular clouds 
(e.g. |Falgarone et al.|2004] ) and may therefore regulate the frag- 
mentation of clouds from large scales down to collapsing cores 
(see Padoan & Nordlund 2002; Cha brier 2003) . The turbulent 
fragmentation can adequately reproduce the observed co re mass 
functions which mimic the IMF of stars ( |Klessen|2001| ). In this 
description, the most massive stars would originate from the rare 
high-mass fragments which could decouple from their turbulent 
environment. This however implies that massive stars would be 
statistically dispersed with large distances between them, in con- 
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tradiction with their observed tendency to cluster. In a purely 
gravitational view, the typical stellar masses are dictated by the 
Jeans masses which are usually of the order of a fraction of a 
solar mass. This can explain the average stellar masses, but has 
problems with explaining the highest masses. The latter are of- 
ten two orders of magnitude larger than the Jeans masses. Some 
additional effects are therefore needed to fully understand how 
fragmentation proceeds to produce massive cores, or to enable 
specific conditions to grow massive stars from low-mass collaps- 
ing fragments. 

The final fragmentation leading to individual collapsing pro- 
tostars inside the dense massive clumps is therefore the central 
question for understanding the origin of massive stars. In the 



slow evolution of turbulent clouds promoted by [McKee & Tan 



( 2002 ) the level of fragmentation is reduced and should lead to 
stable (relatively long living) massive cores which would col- 
lapse into a single star. In the dynamical view, the clumps are 
expe cted to fragment into a large number of low-mass fragments 
(e.g. Bonnell et al.|2007| ) which may grow in mass in a runaway 



process, in the center of the forming cluster when the gravita- 
tional dwell is dominated by the protostellar and stellar compo- 
nents. The fragmentation inside the most massive dense cores at 
the scale of a few 0.1 pc, i.e. the best candidates to be proto- 
clusters in their earliest phases, is highly discriminating for the 
theoretical schemes in competition. 

Since massive stars are rare, and the lifetime of their earliest 
phases is short, the protostellar precursors of massive stars are 
elusive objects. Large scale and systematic surveys are manda- 
tory to find them. The early works obtained in the vicinity of 
luminous IRAS or maser sources in the whole galactic plane 



(Molinari et al.||1996 


|Plume et al. 


1997 


Beuther et al.||2002b| 


Mueller et al. 


2002 


Faundez et a 


l.| 2004 


J), and more recently 



inside the dark filaments seen in the infrared surveys of the 
Galactic plane, the so-called IRDCs (e.g. ^Simon et a l.||2006[ 
and references therei n), have found a number of candidates to 
host such precursors ( Rathborne et al.|2Q06 Pillai et al. 2006). 



These objects are however spread over the whole Galaxy, at dis- 
tances usually larger than 3 kpc, and their detection are not based 
on complete imaging of star- forming complexes, not provid- 
ing a representative view of these not well studied cold phases. 
In Cygnus X, Motte et al. ( 2007| ) could achieve such a com- 



plete survey of a molecular complex rich enough to provide a 
sample of massive dense cores (hereafter MDCs). Cygnus X 
is a nearby mas sive comp lex with a distance of only 1.7 kpc 
dSchneider efal . 2006, S chneider et al.|[2Q07l ). A total of 40 
MDCs has been discovered. These cores are dense, and mas- 
sive on a scale of ~ 0.1 pc and can thus be expected to be rep- 
resentative of the earliest phases of the formation of high-mass 
stars. Even the MDCs which are not seen as dark silhouettes 
on the IR backg round (IRDCs) are re cognized. Following the 
nomenclature by [Williams et al. ( 2000| ), we note that, in terms of 



sizes, the ~ 0.1 pc MDCs of Cygnus X should be considered as 
cores. We here adopt the older and even more classical term of 
"dense cores" for the ~ 0. 1 pc entities even if they might be sub- 
fragmented and may form a small cluster. The term "fragments" 
is then used for smaller structures at the scale of a few 0.01 pc 
inside these dense cores, and which could correspond to individ- 
ual prestellar or protostellar envelopes such as the ones found in 



the ne arby and low-mass clusters in formation (e.g. [Motte et al. 
[1998]). 

We present the first results of a systematic study at high spa- 
tial resolution of the MDCs of Cygnus X. The 1 and 3 mm con- 
tinuum maps of 6 IR-quiet MDCs obtained with the Plateau de 
Bure interferometer are discussed here. The accompanying line 



emission obtained in the same time will be included in forth- 
coming papers dedicated to outflows for the CO and SiO lines, 
and to the dense molecular gas inside the MDCs for the H^^CO^ 
and H^^CN lines. The source selection and the observations are 
explained in Sect. [2] The continuum maps as well as pointed ob- 
servations in N2H^owards each MDCs obtained with the IRAM 
30 m telescope are shown in Sect. [3] In Sect. |4] we discuss the 
obtained properties of the fragmentation of these MDCs in the 
framework of understanding the origin of high-mass stars. 



2. Source selection and observations 

2.1. IR-quiet MDCs in Cygnus X 

From the complete imaging survey of the high column density 
regions of Cygnus X in the 1.3 mm continu um wi th MAMBO 
on the IRAM 30m telescope, Motte_et_alJ p007| ) recognized 



129 dense cores with masses ranging from 4 to 950 M©. 
Among these, a total of 40 were found to be more massive than 
40 M© and were identified to be probable precursors of high- 
mass stars in Cygnus X. With an average size of 0.13 pc, and 
an average density of 1.9 x 10^ cm"^, these MDCs have sizes 
only slightly larger than low-mass dense cores in nearby molec- 



ular clouds ([Ward-Thompson et al.|1994 


Bacmann et al.|2000 


'Evans et al.|2001| [Ward-Thompson et al. 


2007*). But they are on 



average almost 10 times more dense and 20 times more mas- 



sive than their low-mass counterparts (see Tab. 4 in [Motte et al. 
[2007) . Despite such high densities and masses, which should 
lead to intense star formation activity, as many as 17 out of 40 
cores do not emit or emit only weakly in the infrared (IR-quiet). 
This lack of strong infrared emission could indicate that these 
regions were not yet forming massive stars. On the other hand, 
strong SiO emission is systematically detected towards them. 
The intensities and profiles of these SiO lines indicate the pres- 
ence of powerful outflows most probably driven by high-mass 
protostars. These properties are very similar to the properties 
of low-mass Class young stellar objects (hereafter YSOs) as 



dense, cold objects driving powerful outflows (e.g. [Andre et aL 
[1993|). A critical comparison suggests that the IR-quiet MDCs 
are either a scaled-up version (in size and mass) of low to in- 
termediate mass Class YSOs or correspond to clusters with 
a significant number of low-mass Class YSOs. These clus- 
tered Class Os could further compete for mass in the center of the 
cluster if their crossing time is much shorter than their collaps- 
ing time. Alternatively the MDCs may be sub-fragmented into 
a fewer number of high-mass Class YSOs to be identified at 
smaller scales. In all scenarios, these IR-quiet MDCs constitute 
a sample of prime targets for investigating the earliest phases of 
massive star and cluster formation. 

2.2. A representative sample of 6 IR-quiet MDCs 

We have selected the 6 most massive IR-quiet MDCs of Cygnus 
X. The names and coordinates from Motte et al.| ( |2007| ) are listed 
in Tab.JT] An overview of the locations of these MDCs is given 
in Fig.[l[Three of the cores (CygX-N40, CygX-N48, and CygX- 
N53) are located in the DR21 filament which is the most massive 
and dense region of Cygnus X, and a well-known region of mas- 
sive star formation. CygX-N40, and CygX-N48 correspond to 
the submillimeter sources DR21(OH )-N2, DR21( 0H)-S, while 
CygX-N53 is situated close to FIR3 (Chandler et al. 1993 ). The 
other 3 selected MDCs have been discovered by the MAMBO 
survey and are situated in more isolated, but still prominent 
molecular clumps bright in CS emission (Fig. [T]). CygX-N3 in 
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Table 1. Observing parameters and obtained beam and rms. 





Phase center 


Synthesized beam 


Obtained rms 


Source 


(J2000) 


[arcsec] 


[mJy/beam] 


name 


R.A. 


Decl. 


1.3 mm 


3.5 mm 


1.3 mm 


3.5 mm 


CygX-N3 


20 35 34.1 


42 20 05.0 


1.09x0.92 


2.66x2.20 


1.18 


0.212 


CygX-N12 


20 36 57.4 


42 11 27.5 


1.14x0.83 


3.02x2.28 


1.95 


0.221 


CygX-N40 


20 38 59.8 


42 23 42.0 


1.08x0.92 


2.62x2.17 


1.04 


0.171 


CygX-N48 


20 39 01.5 


42 22 04.0 


1.12x0.91 


2.81x2.33 


2.25 


0.340 


CygX-N53 


20 39 03.1 


42 25 50.0 


1.05x0.87 


3.16x2.57 


1.94 


0.246 


CygX-N63 


20 40 05.2 


41 32 12.0 


1.15x0.83 


3.05x2.24 


4.25 


0.332 



FCRAO CS 2-1 

• beam FWHM 




" CygX-N63 



Fig. 1. Overview of the northern part of Cygnus X indicating the 
location of the 6 MDCs from |Motte et al.| ( |20Q7j ). The contour 
map displays the integrated emission in CS(2 - 1) as observed 
with the FCRAO (Schneider et al. 2009, in prep). Note that 3 
of the 6 MDCs (CygX-N40, CygX-N48, and CygX-N53) are 
clustered inside the most massive region of Cygnus X, the DR21 
filament. The red stars indicate the centers of the HII regions 
DR17 to DR23 from |Downes & Rinehart1 ( [T966l ). 



2.3. Interferometric observations of the 1 and 3 mm 
continuum 

We used the IRAlvQ Plateau de Bure Interferometer (hereafter 
PdBI) to image at high resolution the dust continuum simultane- 
ously at 3.5 and 1.3 mm. The 1 mm and 3 mm receivers were 
tuned to 230.54 and 86.61 GHz respectively with two 320 MHz 
units placed at 230.41 and 230.67 GHz and 86.48 and 86.74 
GHz respectively. In addition to the continuum, the emission in 4 
molecular lines has also been correlated: ^^C0(2- 1), SiO(2- 1), 
Hi3C0+(l - 0), and H13CN(1 - 0) at 230.54, 86.85, 86.75, and 
86.34 GHz respectively. 

The observations were done in track-sharing mode with two 
targets per track with the following pairs: CygX-N48/CygX- 
N53, CygX-N3/CygX-N40, and CygX-N12/CygX-N63. The D 
configuration tracks were obtained between June and October 
2004 with 5 antennas with baselines ranging from 24 m to 82 
m. The C configuration tracks were obtained in November and 
December 2004 with 6 antennas in 6Cp with baselines ranging 
from 48 m to 229 m. We mostly used as phase calibrator the 
bright nearby quasar 2013-^370 and as flux calibrator the evolved 
star MWC349 which is located in Cygnus X. 

In order to favor the spatial resolution, the maps were 
cleaned using the uniform weighting. The resulting synthesized 
beam and rms in the continuum are summarised in Tab. [T] to- 
gether with the field names and centers of phase. The clean- 
ing components were searched in the whole area of the primary 
beams. No polygon was used in order to avoid introducing any 
bias in the resulting maps of emission. 

2.4. I RAM 30 m observations of N2H^ 1^0 

The N2H^ 1^0 line at 93.176265 GHz was observed in Septem- 
ber 2003 using the BlOO receiver and the VESPA correlator at 
the IRAM 30m telescope. The average system temperature was 
182 K and the average rms of the spectra is 0.19 K. The data 
are corrected for the main beam efficiency of 0.78 and have a 
velocity resolution of 0.13 km s"\ and a beamsize of 29''. 



the western part, is located close to DR17 which is an HII region 
excited by two OB clusters (clusters #12 and 14 in'Le Duigou 
& Knodlseder||2QQ2| ). These clusters are shaping the cloud and 
CygX-N3 corresponds to the tip of a pillar-like cloud. CygX- 
N12 is also situated in a cometary shape cloud probably influ- 
enced by the same OB clusters but the cloud seems to be less 
compressed from outside than CygX-N3. Finally CygX-N63 is 
placed in the south of DR21 in the DR22-DR23 filament (see 
[Schneider et al.|2Q06) . It can be considered as the most isolated 
MDC of the sample. 



3. Results 

3.1. The 3mm continuum maps 

The left panels in Figs.[2]and[3]display the dust continuum emis- 
sion towards the six cores at 3.5 mm. Strong 3 mm emission is 
detected for all fields except CygX-N40. Despite a total mass of 
106 Mq in Motte et al., CygX-N40 seems to be dominated by ex- 
tended emission which is fully filtered out by the interferometer. 
We note that it has the lowest average density in our sample (see 



1 IRAM is supported by INSU/CNRS (France), MPG (Germany) and 
IGN (Spain). 
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Fig. 2. Maps of the continuum at 3.5 mm (left) and 1.3 mm (right) obtained with the PdBI for the three more isolated MDCs (see 
text) . The color images have linear scales from -2cr to +20 cr for the 4 top panels, and from -3 cr to +30 cr for the 2 lower ones 
where cr is the resulting rms in each field (see Tab. [T])(the unit on the right of the plots is Jy/beam). The (logarithmic) contour 
levels are 3cr, 4.8 cr, 7.5 cr, 12 cr, 19 cr, 30 cr, 75 cr, and 120 cr. The crosses and numbers indicate the locations and the names of 
the fragments at 1 mm (see Tab.|3]). The synthetised beams are displayed in the bottom right corners. The dashed circles indicate 
the primary beams. In the 3.5 mm maps (left), the name of the MDC is given, and the dashed squares indicate the areas covered at 
1.3 mm (right). Note the similar behavior of the 3 MDCs. A compact core is observed at 3.5 mm which splits into smaller fragments 
at 1.3 mm. CygX-N63 seems to actually correspond to a single object down to the r'(1700 AU) spatial resolution. 



Tab.|2]), and it is located in the large scale DR21 filament where 
it is possible that the observed emission is more elongated in 
the direction of the line-of-sight. This could have led us to mis- 
interpretate CygX-N40 as a single MDC and to over-estimate the 
local density. It will therefore not be further considered for dis- 
cussion. For the other fields, the emission contains bright com- 
pact emission peaks. They are structured into a few individual 



bright compact sources - except CygX-N63, which seems to be 
dominated by a single object. We note also that CygX-N48 ap- 
pears less centrally concentrated than the other cores, and could 
be dominated by a large population of unresolved compact frag- 
ments forming a full cluster. 
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Fig. 3. Same as Fig. [I for the three MDCs located in the filament of DR21 (see text). Note that in contrast to the isolated MDCs 
(Fig. [2]), the cores show here diverse structures. CygX-N40 appears devoid of compact emission at both 3.5 and 1.3 mm. CygX- 
N48 has a complex and extended structure at 3.5 mm with two main cores in the central regions. The 1.3 mm emission resolve 
a few, probable compact, fragments in the two main cores. The bright peak at 3.5 mm in the South-East corresponds to the weak 
1.3 mm source outside in the primary beam in the right panel. It is located at the edge of the nearby cluster associated with 
MSX6G81. 7133+0. 5589. CygX-N53 has a bright central core at 3.5 mm which splits into 4 smaller fragments at 1.3 mm in a way 
similar to the isolated MDCs. The 3.5 mm source in the SE is hardly detected at 1.3 mm due to the falloff* in sensitivity of the 
primary beam (dashed circle). 



3.2. The 1mm continuum maps 

The right panels in Figs. |2] and [3] display the dust continuum 
maps at 1 mm. For three of the cores (CygX-N3, CygX-N12, 
CygX-N53) the diff'erence between low (3.5 mm) and high 
(1.3 mm) resolution maps is striking. While at a scale of 2''to 
3 ''resolution, the emission mainly consists of extended, elon- 



gated structures, at the ~ T'resolution of the 1.3 mm map the 
emission clearly splits into a number of individual usually spher- 
ical fragments. In contrast, CygX-N48 has a similar aspect than 
at 3.5 mm with emission spread over the whole region and with 
compact sources which are not very spherical such as if it con- 
sists of partly non resolved individual objects. CygX-N63 seems 
to diff'er from the others since the emission is dominated at 
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Fig. 4. N2H^(1 - 0) spectra towards the center of each MDCs observed with the IRAM 30 m. Note that the velocity scale is given 
for the weakest line (left component) which is the only isolated component of the hyperfine structure. The other 5 components of 
the hyperfine structure are blended and form the 2 apparently brighter groups of lines. 



all observed scales by a single central object. A slight trend 
of having some distributions of aligned compact sources inside 
filamentary structures seen at 3.5 mm can be noted. CygX-N3 
and CygX-N12 are composed of three to four aligned compact 
sources. CygX-N48 has a "S" shape emission which could be a 
filamentary structure too. CygX-N53 displays an almost north- 
south ridge of emission which is aligned with the large scale ori- 
entation of the DR21 filament in which it is embedded. The di- 
rection of the axis between the two brightest sources is however 
oriented east- west, almost perpendicular to the main direction of 
the large scale emission. 



3.3. Source extraction at 1 mm 



Our ultimate goal is to study the population of individual pro- 
tostellar objects inside the MDCs of Cygnus X. We therefore 
extract fragments in the highest spatial resolution data maps 
(1 mm maps). To be systematic, we used the Gaussclump algo- 
rithm ( Stutzki & Guesten 1990) adapted to continuum maps as 
described in Motte et al., ( ,2007 ) to extract the compact sources. 
Above a safe 7 sigma level, 24 sources have been detected inside 
the six fields and 23 inside 5 fields (excluding CygX-N40), i.e. 
an average of almost 5 sources per field. All sources are clearly 
identified already by eye-inspection in the maps. 

The obtained list of fragments, ordered by increasing peak 
flux, is given in Tab. [3] The table gives the coordinates, peak 
fluxes, 2D gaussian sizes (FWHM), and integrated fluxes di- 
rectly obtained from the Gausclump analysis. The integrated 
fluxes are converted into masses Mphrt, using a d ust emissivity 
at 1.3 mm of 1.0 cm^g"^ ( Qssenkopf & Henning 1994), a gas-to- 
dust mass ratio of 100, and a dust temperature of 20 K (formula 
(1) in Motte e t al. 200 7). The de-convolved sizes (col. 8 in Tab. [3]) 
correspond to the quadratic average of the de-convolved 2D sizes 
(FWHM) and are converted in AU for d=1.7 kpc. Average densi- 



ties are calculated inside the projected region of radius Rgo con- 
taining 90 % of the mass (MpdBi) for a 2D gaussian distribution: 



Rgo = 1.95 X FWHM / V81n2. 



(1) 



The masses MpdBi obtained from the gaussian decomposition 
of the PdBI maps are good estimates of the masses inside the in- 
ner regions of the compact fragments. Their actual sizes and total 
masses can however not be easily derived from the present inter- 
ferometric data which tend to filter out the extended envelopes 
(see Sect. 4.3 ). An estimate of the total masses Mtot is given in 
Tab. [3] This is an extrapolation of MpdBi inside a total size of 
3000 AU (FWHM). According to the detailed discussion given 
in Sect. |4.3| these masses should be seen as conservative esti- 
mates of the total masses in contrast to MpdBi which are only 
lower limits. 

We take a temperature of 20 K for all mass derivations fol- 
lowing Motte et al. ( |2007| ). This typical temperature has b een 



recently confirmed by NH3 measurements by Wienen|p008 



ing the Eff'elsberg t elescope. The emissivity adopted (Ossenkopf] 



us- 



& Henning] 1 1994| ) is used in most studies of low-mas s star- 
formation but also for high-mass protostellar object s ([Andre 



eF^"2000 , H unter et al.||I5Q6| [Rathborne et"aL||2507| [SchuUer 
et al.|,200 9) and should be adequate for dense an d cold gas. 
Alternatively, some authors use the older values by [Hildebrand 



(1983) usually in the case of warm objects for which the emis- 
sivities might be slightly difl'erent. These emissivities are 3 times 
smaller than the ones of ||Ossenkopf & Henning ( 1994 ) in the mil- 
limeter range which is a significant diff'erence, leading to masses 
~ 3 times larger than the ones we derived here. Any compar- 
ative analysis should therefore take into account these diff'erent 
assumptions on emissivities. 

3.4. 1^0 emission towards the MDCs 

The obtained single spectra towards each MDC are shown in 
Fig.|4] We use these spectra to estimate the average turbulent ve- 



S. Bontemps, R Motte, T. Csengeri, N. Schneider: Fragmentation in the massive dense cores of Cygnus X 



Table 2. Global properties of the MDCs derived from single dish 
continuum and N2H^ line observations. 



Name 


Mass 


Size 


< «H2 > 


Cturb 




■^Jeans 


-^Jeans 


^MDC 


Tff 






[Mo] 


[AU] 


[cm~^] 


[km.s-1] 


[Mo] 


[AU] 


[Mo] 


[g.cm-2] 


[yr] 


[AU] 


CygX-N3 


84.0 


20300 


5.29 X 10' 


0.98 


91.0 


8490 


1.1 


1.45 


6.0 X 10^ 


6150 


CygX-N12 


86.0 


20400 


5.33 X 10^ 


0.89 


75.5 


8460 


1.1 


1.47 


5.9 X 10^ 


3420 


CygX-N4(f] 


106.0 


33200 


1.52 X 10^ 


0.82 














CygX-N4r 


197.0 


27300 


5.10 X 10^ 


1.28 


209 


8640 


1.1 


1.88 


6.1 X 10^ 


3830 


CygX-N53 


85.0 


24000 


3.24 X 10^ 


0.76 


64.8 


10850 


1.4 


1.05 


7.6 X 10^ 


4180 


CygX-N63 


58.0 


1230(f] 


1.64 X 10^ 


0.72 


29.8 


4820 


0.62 


2.73 


3.4 X 10^ 





^ CygX-N40 does not show any bright compact emission in the PdBI 
maps suggesting it has been misinterpreted as a MDC probably due to 
confusion on the line-of-sight towards the DR21 filament (see text for 

details). It is therefore not further consider for discussion. 

^ CygX-N63 was actually not resolved out by MAMBO in 



Motte 



|et al.| ( [2007] ), and this size might actually correspond to an upper limit 
of the source size. 



locity dispersion in the cores. N2H^ 1^0 is well suited since it 
is usually optically thin and because its abundance drops dras- 
tically in the lower density, warmer outer layers which tend to 
dominate the line widths of most other molecular tracers. Since 
the profiles are not always perfectly gaussian, we did not use 
gaussian fits, but measured directly the widths at half maximum 
of the lines. For these measurements, we only use the isolated 
component of the hyperfine structure (HFS) which appears as 
the weakest line in Fig. [4] The other 5 components are blended 
and form the 2 brighter groups of lines. The resulting values con- 
verted into the gaussian ID velocity dispersion (Xturb are given in 
Tab. [21 

The fits of the 6 components of the N2H^ 1^0 multiplet 
can provide some indications of the opacity and of the "aver- 
age" excitation temperature. The 6 relative distances and inten- 
sities of the HFS components of the 1^0 line are given as in- 
put to a simultaneous gaussian fit to all components, assuming 
equal excitation temperature, using the 'Method HFS" feature of 
GILDAS|5 These fits confirm that the lines are actually optically 
thin with values of r between 0.03 and 0.2. Also the obtained 
excitation temperatures are typically around 20 K, conform with 
NH3 observations. 



inside R90 (Equ.[T]), the radius which corresponds to 90 % of the 
mass for the 2D gaussian distribution used to extract the cores. 

In the second part of Tab. |2] we give the quantities which 
are believed to ultimately dictate the properties of the proto- 
stellar fragmentation and of the subsequent star formation: the 
virial masses Myir, the local Jeans lengths /Ijeans and Jeans masses 
^jeans, the mean surface densities Smdc, and the free-fall times 
Tff. In addition the typical, observed separations /Isep between 



the PdBI fragments is provided in the last column (see Sect. 4.4 
for details). The Jeans masses, Jeans lengths, and virial masses 
are calculated for the average conditions in the MDCs using the 
classical formulation: 



'i^Jeans 



Jeans 




47r (Aj 
-Po 



(2) 



(3) 



(4) 



4. Discussion 

4.1. Global properties of the MDCs 

The targeted MDCs are very dense cores with densities 10 times 
larger than i n low-mass star-forming regions (see Tab. 4 in Motte] 
et al.||20Q7] ). They have total masses large enough to be single 
high-mass protostars or to form a whole cluster of stars. We re- 
view here their global properties in order to discuss which popu- 
lation of stars these cores are expected to form, based on the the- 
oretical predictions of the usually promoted scenarios for high- 
mass star formation. 

In the first columns of Tab. [2] we give the properties of the 
MDCs as directly derived from the MAMBO survey of Motte 
^taL] ( |2QQ7l ): the names, total masses, ID deconvolved sizes, and 
average densities, as well as the ID turbulent dispersion derived 
from the N2H^ line widths. The average densities are however 
not exactly the same as in |Motte et al. ('2007 ), since we derive 
them in slighly diflTerent volumes. We average here the densities 



which numerically gives 



..a. = 1.38xlO^Au( — ) (^) 



Grenoble Imaging and Line Data Analysis software 



M... = 0.90Mo(4)"(^) 



-0.5 



115.6M,r 



R 



0.1pc/\lkm.s 



<^turb 



(5) 
(6) 
(7) 



The listed /Ijeans and Mjeans have been calculated for a tempera- 
ture of 20 K (see Sect. |3.3| ). In an earlier phase for the real initial 
conditions, the MDCs could have been colder, perhaps as down 
as 10 K such as what is observed in low-mass dense cores, lead- 
ing to smaller Jeans lengths and masses. The virial masses are 
calculated for R = Rgo and (Xturb is the ID (average on the line- 
of-sight) velocity dispersion. The values for Smdc = M/(nR^) 
are calculated for the h alf-mass radius for pr oper comparisons 
with values promoted in McKee & Tan ( 2003| ). 
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The average densities in the MDCs are so high that the free- 
fall times are of the order of or smaller than 10^ yr. This is ac- 
tually smaller than the statistical timescales for low-mass proto- 
stars derived in nearby star-forming regions. The duration of the 
protostellar evolution, which corresponds to the sum of the life- 
time s of Class and Class I YSQs, is of the order of 2 x 10^ yr 
(e.g. Greene et al ."1994i Ken yon & Har tmann'1995 ). The free- 
fall times are also roughly equal to the crossing times of the 
cores, estimated by dividing the sizes by the 3D velocity disper- 
sions ( V3 X (Tturb). which are typically equal to 6 x 10"^ yr. The 
MDCs are clearly self-gravitating since the ratios M/Myir are all 
close to or even larger than 1 and up to 2 for CygX-N63. The tur- 
bulent support can therefore not prevent the cores from collapse. 
The unknown magnetic field could delay the collapse but cannot 
avoid it. If they actually collapse in a free-fall time, the global in- 
fall rates will be of the order of a few 10"^ Mo/yr which is high 
enough to overcome the radiation pressure barrier. These MDCs 
could therefore be in monolithic collapses leading to single mas- 
sive stars. Their global properties (total masses, sizes, turbulent 
velocity dispersion, and mean surface densities) are conform 
with the initial conditions for the turbulence-regulated mono- 
lithic collapse proposed by M cKee & Tan| ( |2002| ). Alternatively, 
according to |Dobbs et al.^ ( 2Q05) for instance, such massive, tur- 
bulent cores should actually fragment to form a whole cluster of 
stars. 



4.2. Properties of low-mass, nearby protoclusters 

The earliest phases of low-mass star formation as observed in 
nearby star-forming regions correspond to the so-called pre- 
stellar cores and to Class and Class I protostars. They ap- 
pear to form in clusters in regions such as p Ophiuchi, Perseus, 
Serpens, or Orion. Below we refer mostly to results obtained in 
p Ophiuchi and Orion because these are reference regions, and 
our group has studied them using the same, homogeneous way 
to derive masses, and sizes like in Cygnus X. References for 
other nearby protoclusters can be found in recent r e views such as 
Andre et al. ( |2QQQ ), Ward-Thompson et al. ( 2007| ), di Francesco 



etal.| ( |2Q07| ) 



One of the important findings obtained in these regions is 
that the fragmentation of clouds and dense cores is probably at 
the origin of the stellar masses, and that the masses of objects at 
the beginning of the collapse accurately reflect the final stellar 
mass in absolute values, or at least recognizing that a large frac- 
tion (efficiency) end up in the star. Motte et al. ( 1998 ), followed 
by a large number of other works (e.g. |Tes ti & Sargent 1998( 
Johnstone et"ar|[2000[ [Nutter & Ward-Thompson 2007) , could 



trace that the mass distribution of pre-stellar cores and Class 
YSOs have the same shape than the IMF. The stellar masses 
seem to be determined already at the pre-stellar stage as a re- 
sult of the fragmentation of dense cores. Inside these low-mass 
dense cores, the fragmentation eflftciency seems to vary with the 
average density and is typically of the order of 5 to 20% but has 
been found to reach 40% in p Oph A (see Fig. 7 in Motte et al. 
[1998 ). The fragmentation efficiency is defined as the fraction of 
mass inside the pre- and proto-stellar envelopes. On the other 
hand, the global star formation efficiency (hereafter SFE), which 
includes also the masses of the already formed stars , is on aver- 
age more of the order of 30% ( Bontem ps et al.|20QT] ) at the scale 
of the dense cores in the p Ophiuchi proto-cluster. The typical 
sizes of individual protostellar objects are of the order of a few 
1000 AU (3000 AU in p Ophiuchi and 5000 AU in Orion B for 



instance) and the separations between objects was found to be 
typically of 5000 AU in nearby proto-clusters. 

4.3. A population of dense, self-gravitating fragments 

The MDCs in Cygnus X are clearly massive and self-gravitating. 
They are therefore expected to contain collapsing objects which 
could be true massive protostars. Our new PdBI maps at high 
resolution down to physical scales of only 1700 AU should re- 
veal these protostellar objects. 

At the moderate spatial resolution of ~ 3'' (5 100 AU) at 
3.5 mm (left panels in Fig. [2] and |3]), 4 out of the 5 MDCs con- 
sist of strong, centrally concentrated, emission emanating from 
single bright, sometimes elongated (CygX-N3 and CygX-N12), 
cores in the inner regions of the MDCs. But, at the highest spatial 
resolution of ~ r'(1700 AU) at 1.3 mm (right panels in Fig.|2] 
and[3]), the picture drastically changes. The inner cores of CygX- 
N3, CygX-N12, and CygX-N53 actually split into a number of 
roughly spherical fragments. A total of 11 of such fragments (3 
to 4 per MDC) are found to have roughly similar sizes and sep- 
arations. Only CygX-N63 stays singlq^in the 1.3 mm map. The 
case of CygX-N48 appears slightly diff'erent with emission at 
3.5 mm which is significantly less centrally condensed than the 
other 4 MDCs. The central emission already splits at 3.5 mm 
into at least two main cores instead of one in the other MDCs. 
At 1.3 mm, 4 possibly individual fragments could be resolved 
and extracted in two main cores, but these fragments stay less 
spherical and with smaller separations than in the other MDCs. 
The diff'erence may be due to a less centrally condensed MDC, 
and to a possible larger population of weaker fragments which 
are not fully resolved due to crowding. Finally, CygX-N63 is sin- 
gle and thus counts for one additional fragment. Altogether, we 
therefore consider that a population of 16 individual fragments 
has been recognized in the central regions of 5 MDCs. The re- 
maining 7 sources listed in Tab 3] are weaker sources usually in 
the outskirts of the MDCs which can be considered as more un- 
certain detections (some are not confirmed at 3.5 mm). 

We see that 2 out of the 23 individual fragments have MpdBi 
larger than 8 M0(CygX-N53 MMl and CygX-N63 MMl), and 
5 additional ones larger than 4 M©. The values MpdBi actually 
measure more the masses of the inner regions of the fragments 
than their total masses (see discussion below). Such masses are 
among the largest masses ever measured in molecular clouds, 
inside physical sizes of only -1000 AU. The resulting average 
densities are extremely high, up to more than 10^ cm"^ (mea- 
sured in Rgo) for 5 of them, and the resulting free-fall times are 
small, with an average value of only 3900 yr for the inner re- 
gions of these 7 fragments. The typical velocity dispersion in the 
inner regions of the MDCs is not known. It is however reason- 
able to consider that it is not larger than the average dispersions 
cTturb measured at the scale of the MDCs. In fact, the velocity 
dispersion is expected to decrease from large to small scales. If 
it would follow a classical Larson's law of index 0.5, the ve- 
locity dispersion, which is on average equal to ~ 0.9km.s"^ at 
the 0.1 pc scale, would be 4.5 smaller at the 1000 AU scale, 
i.e. of the order of 0.2km.s"^ which is close to the thermal ve- 



^ The secondary source CygX-N63 MM2 does not seem to be con- 
vincingly separated from MMl. The comparison of the 3.5 and 1.3 mm 
maps shows that the global shape of the emission is unchanged with 
clear extensions of MMl towards the south-west, east and north-west 
directions at both resolutions and wavelengths. The extraction of MM2 
as an individual fragment seems more due to this east extension than to 
the existence of a secondary object. 



S. Bontemps, R Motte, T. Csengeri, N. Schneider: Fragmentation in the massive dense cores of Cygnus X 9 



Table 3. List of fragments and physical properties at 1.3 mm. 



Fragment 


Coordinates 


Peak 


FWHIVJ^ 


Int. 


MpdBI 


Deconv 


PI — 


A yr 1/71 

MenvQl 


Name 


R.A. 


Decl. 


Flux 


Flux 




Size 






(J2000) 


(J2000) 


[mJy/beam] 




[mJy] 


[Mo] 


[AU] 


r -In 

[cm ^] 


[Mo] 


CygA-JNJ MMl 


20 35 34.63 


42 20 08.8 


O'? o 


1 T7 1 rvc 


/in ^ 


O /CO 

Z.DJ 




1 /I o 1 r\8 
1.4o X lU 


o.Ud 


CygA— iNJ MMz 


20 35 34.41 


42 20 07.0 


dKj.Z 


i.J4 X i.U4 


/1 1 ^ 


z.zi 




1 ni \/ 1 n8 
i.Ui X iU 


^ Qi 


CygX-N3 MM3 


20 35 34.23 


42 20 04.7 


14.4 


1.01 X 1.00 


14.5 


0.77 


< 851 


>6.56 X 10 


2.71 


AT" ATO A TA It A 

CygX-N3 MM4 


20 35 34.55 


42 20 00.3 


8.7 


2.01 X 1.01 


17.7 


0.94 


1720 


9.65 X 10" 


1.63 


CygX-N12 MMl 


20 36 57.65 


42 11 30.2 


64.1 


1.22 X 1.12 


92.6 


4.91 


1100 


1.93 X 10^ 


13.36 


CygX-N12 MM2 


20 36 57.51 


42 11 31.2 


54.7 


1.46 X 1.13 


95.2 


5.05 


1420 


9.23 X 10^ 


10.65 


CygA-iN iz MM J 


20 36 57.81 


42 11 29.7 


10/1 

iz.4 


1 in \/ 1 1 Q 

i. / / X L.Ly 


Z / A 


1 /I ^ 

1.4 J 


1 Qon 


1 0^ \/ 1 n7 

i.ZD X iU 


Z.jy 


/^»,~AT" ATI'^ A TA T /I 

CygX-N12 MM4 


20 36 57.45 


42 11 32.4 


11.7 


3.11 X 0.98 


37.6 


1.99 


2460 


7.00 X 10 


2.42 


CygX-N40 MMl 


20 38 59.54 


42 23 43.6 


5.7 


1.04 X 1.00 


5.9 


0.31 


< 847 


>2.72 X 10^ 


1.11 


CygX-N48 MMl 


20 39 01.34 


42 22 04.9 


53.2 


1.58 X 1.07 


88.5 


4.69 


1400 


9.05 X 10^ 


10.07 


CygX-N48 MM2 


20 39 01.10 


42 22 08.3 


37.5 


1.84 X 1.28 


86.6 


4.59 


1960 


3.19 X 10^ 


7.01 


CygX-N48 MM3 


20 39 01.31 


42 22 02.8 


25.6 


1.65 X 1.01 


41.7 


2.21 


1360 


4.63 X 10^ 


4.88 


CygX-N48 MM4 


20 39 01.25 


42 22 08.3 


20.0 


1.64 X 1.14 


36.8 


1.95 


1570 


2.65 X 10^ 


3.72 


CygA-JN45 MM J 


20 39 01.92 


42 22 02.9 


1 o o 

iz.o 


O O /I 1 O /I 

Z.j4 X 1.Z4 




1 no 


OOO A 


<5.U4 X lU 


O /I o 

Z.4o 


CygX-N53 MMl 


20 39 02.96 


42 25 51.0 


106.3 


1.19 X 1.09 


151.0 


8.00 


1050 


3.61 X 10^ 


22.80 


CygX-N53 MM2 


20 39 03.22 


42 25 51.2 


71.1 


1.39 X 1.03 


112.1 


5.94 


1230 


1.66 X 10^ 


14.44 


L.ygA-iNjJ MM J 


20 39 03.05 


42 25 49.9 


1 T 


i.4D X 1.4 J 


/in 1 


O 1 


1 cm 
io /U 


i. /z X iU 


'J /lO 

J.4z 


CygX-N53 MM4 


20 39 03.02 


42 25 52.7 


16.2 


1.44 X 1.44 


37.0 


1.96 


1840 


1.66 X 10^ 


3.20 


CygX-N53 MM5 


20 39 02.76 


42 25 43.3 


13.2 


1.72 X 1.60 


39.8 


2.11 


2300 


9.08 X 10^ 


2.75 


CygX-N53 MM6 


20 39 02.93 


42 25 45.3 


11.8 


2.21 X 1.02 


29.4 


1.56 


1980 


1.06 X 10^ 


2.36 


CygX-N53 MM7 


20 39 02.87 


42 25 46.9 


9.5 


2.27 X 1.07 


25.3 


1.34 


2100 


7.60 X 10^ 


1.92 


CygX-N63 MMl 


20 40 05.39 


41 32 13.1 


274.3 


1.39 X 1.05 


419.2 


22.22 


1210 


6.64 X 10^ 


55.20 


CygX-N63 MM2 


20 40 05.51 


41 32 12.7 


44.1 


2.09 X 1.44 


139.4 


7.39 


2440 


2.68 X 10^ 


9.09 


CygX-N63 MM3 


20 40 04.97 


41 32 22.5 


37.8 


1.63 X 1.49 


96.1 


5.09 


2060 


3.05 X 10^ 


7.41 



^ Not deconvolved FWHM of the 2D gaussian. 

^ Deconvolved FWHM in ID (quadratic average of the 2 axis sizes) 
expressed in AU for d=1.7 kpc. 

The average density is calculated inside the projected region con- 
taining 90 % of the mass in a 2D gaussian assuming a spherical distri- 
bution. 

MpdBi extrapolated from the deconvolved size to a FWHM size of 
3000 AU assuming a density profile in r~^. 



locity dispersion (0.27 km.s"^ for 20 K). Even in the extreme 
case there would be no decrease at all of the velocity dispersion 
at small scales, with a velocity dispersion of 0.9km.s"^ in the 
inner regions, the virial mass for Rgo = 830 AU (which corre- 
sponds to a size of 1000 AU in Tab|3]) is 3.8 Mq which is still 
smaller than IMpdBi for these 7 most massive fragments. We can 
therefore safely conclude that these fragments, or at least their 
inner regions, are self-gravitating, and are most probably col- 
lapsing objects. Another indication of the probable collapsing 
nature of these fragments is the general trend of the fragments 
to be spherical in the continuum maps such as if gravitation is 
shaping them. 

The typical separations between the fragments can be evalu- 
ated from the projected distances between the fragments. The de- 
projectec^ average separations /Isep between the closest neigh- 
bors of the 16 main fragments are given as the last column of 
Tab. [2] We note that the separations are of the order of 5000 AU, 
typical for values derived in nearby, low-mass proto-clusters. 

We have thus detected a population of self-gravitating frag- 
ments whose sizes and separations are similar to low-mass pre- 
stellar cores and protostars in nearby proto-clusters. These frag- 
ments are most probably individual protostellar objects of sim- 
ilar nature, but more massive, than in nearby regions. Since the 



^ We use a de-projection factor equal to l/sin(57.3°)=1.19 which cor- 
rects for the unknown random orientation of the separation axis. 



selected ]V[DCs stayed undetected in the infrared, most of these 
objects should still be cold and young and should thus be more 
pre-stellar cores or Class YSOs than Class I YSOs. Also the 
selected IMDCs are not associated with any strong radio centime- 
ter sources from the literature. Three of them, CygX-N3, CygX- 
N12, and CygX-N63, were confirmed to be devoid of 3.6 cm 
continuum sources at low levels, down to 0.1 mJy thanks to re- 
cent dedicated VLA observations (Bontemps et al. in prep). No 
UCHii region is therefore yet expected inside these cores. Also, 
several collimated CO outflows have been found to be driven 
by most of the massive fragments in our accompanying C0(2-l) 
PdBI data. These data will be presented in a forthcoming paper 
entirely dedicated to the nature and evolutionary status of the 
protostars. 

4.4. Total masses of the protostellar fragments 

The sizes and masses of the fragments derived from the interfer- 
ometer maps are aflTected by the spatial filtering due to the lack of 
information for the large scales (typically larger than 1/3 of the 
primary beam) but also due to the limitation of the deconvolu- 
tion algorithm for limited UV coverages and noisy phases at high 
frequencies (1.3 mm) which aflfects also the intermediate spatial 
scales (sidelobes which cannot be cleaned). The stronger spatial 
filtering at 1.3 mm compared to 3.5 mm actually partly explains 
why the individual small scale fragments show up so clearly at 
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1.3 mm, the different spatial resolution does not explain all the 
differences. It is particulary visible for CygX-N3 and CygX-N53 
where the 3 mm maps show strong, central, slightly extended 
emissions which do not seem to be entirely recovered at 1.3 mm 
where only the compact fragments are well detected. We there- 
fore expect that the masses obtained at 1.3 mm are only valid for 
the smallest scales. We should thus regard the masses MpdBi and 
sizes of the fragments as lower limits. 

In order to estimate what could be the true total masses of 
the fragments, we have extrapolated MpdBi up to realistic sizes 
for such proto-stellar objects, and assuming a classical density 
profile in r"^. The most reliable sizes are provided by what has 
been observed in nearby protoclusters. In the p Ophiuchi proto- 
cluster, [Motte et al.| ( |1998| ) derived an average size for 58 starless 
millimeter condensations of the order of 3000 AU (deconvolved 
FWHM) for masses in the range 0.1 to 1 M©. In Orion B, the 
more massive condensations (masses up to ~3 M©) were found 
to be slightly larger with an average size closer to 5000 AU (see 
Fig. 4 in Motte et al.|2001| ). By comparing p Ophiuchi and Orion, 
it seems that the cores tend to be larger in size when they are 
more massive. In Cygnus X, and for higher mass pre-stellar and 
protostellar cores, we could therefore expect even larg er sizes for 
the collaps ing cores. This is also what is predicted in |McKee"& 
|Tan| ( |2003| ). We therefore adopt here the conservative, smaller 
size of 3000 AU like in p Ophiuchus to extrapolate the PdB 
masses, and to derive the values Mtot given in Tab. [3] 

4.5. Progenitors of OB stars ? 

The obtained masses Mtot then represent the gravitationally 
bound masses associated with each proto-stellar object. The bulk 
of these masses is expected to end up in the final star. Only the 
feedback from the outflow can decrease the global efficiency. 
Since outflows are well collimated during the earliest phases of 
the accretion phases (Class Os) that corresponds to the phase for 
which most of the mass is accreted (e.g. [Bontemps et al.|1996 ) 
it is reasonable to assume that the efficiency stays rather high 
even considering outflows. The outflows from high-mass Class 
Y SOs are less known and could perhaps be less collimated 
(e.g. Arce et al. 2007 ). To be cons ervative, we thus as sume here 
a rather low efficiency of 70 % (e.g. [Bontemps et al.|2001| ). Using 
this efficiency, we derive that all fragments with initial mass 
larger than 11.4 MqwIII form a star more massive than 8 M©. 
This is the case for 4 fragments. But, of course, statistically, we 
expect that the observed protostars have already accreted a frac- 
tion of their final mass. With a typical fraction of 50 % of the 
mass already accreted, the statistical limit between intermediate 
mass protostars and high-mass protostars is then equal to 5.7 M©. 
This is the case for 5 additional fragments. It is therefore a total 
of 9 fragments which could be precursors of ionizing OB stars 
(earlier than B3, > 8 M©). This estimate is based on statisti- 
cal arguments and is therefore subject to uncertainties. The total 
number of O stars in the whole Cygnus X region is necessary 
larger than the 100 O stars recognized in the 2 Myr old Cygnus 
0B2 association (see Knodlseder[[2000[ ). Assuming 50 % more 
O stars spread over the whole Cygnus X complex, we expect 
150 O stars in total, and of the order of 750 ionizing OB stars 
(assuming a Kroupa 2001 IMF) formed during the last 2 Myr. 
For a constant star formation rate, the 9 OB star precursors in 
the MDCs of Cygnus X would then correspo nd to a timescale o f 

2.4 X 10^ yr quite typical for Class YSOs ( Andre et al. 2000). 
This time is compatible with the free-fall times of Tab.[2]and cor- 
responds to 4 Mq (50 % of the lowest stellar mass of OB stars) 
accumulated at a rate of 1.7 x 10"^ Mo/yr. 




LOG(M) 



Fig. 5. Mass function of the 23 fragments in the 5 MDCs consid- 
ered in the discussion (grey histogram) using the values Mtot- For 
comparison the global stellar mass function is shown for a global 
star formation efficiency of 30 % of the total mass of 51OM0for 
the 5 MDCs, and using the Kroupa (2001) IMF for between 
0.01 and 12OM0(blue histogram in the background). The low- 
mass fragments in the MDCs cannot be detected but it seems 
that already the detected massive objects are over-massive com- 
pared to what a normal IMF would predict. Even in the extreme 
view of using the highly filtered masses MpdBi, there is still a 
slight excess of massive objects (dashed blue histogram). 



4.6. Comparison with previous studies of MDCs 

A number of high-mass star-forming regions has already been 
investigated at high spatial resolution with (sub)millimeter inter- 
ferometers, mostly with the IRAM PdBI and the SMA. But only 
a few studies were dedicated to the earliest phases of star for- 
mation and to regions of 0.1 pc or smaller. In fact, the presently 
available millimeter interferometers are limited in spatial reso- 
lution and only the most nearby MDCs can provide access to the 
smallest physical scales required to reach the individual proto- 
stellar objects. A few nearby HMPOs from the Sridharan et al. 
(200 2| sample were o bserved with the PdBI (see for instance 
Beuther et al.[[2002at [Leurini et al. 2007), while a sam ple of 
not too distant IRDCs is presented in Rathborne et al.[ ( |2007[ ). 
These MDCs were all found to be sub-fragmented, with the pres- 
ence of a few high-mass fragments or condensations which are 
good candidates to be direct precursors of OB stars. A mass 
fu nction of intermediat e mass fragments could even be derived 
by [Beuther & Schilke[ ( [2004] ) for IRAS19410-h2336, where it 
is recognized that the core mass function has roughly an IMF- 
like shape towards the high-mass regime too. However, the de- 
tected high-mass condensations were not so clearly recognized 
as compact, spherical fragments. Detections show only elon- 
gated emission distributions and weak peaks, moreover, many 
at the limit of the detection level (e.g. the PdBI 1 mm maps 
of [Rathborn e et al.|[200 7). The better q uality of our m aps com- 
pares more with the results obtained by Peretto etliL] ([2007 ) for 
a nearby MDC (d = 800 pc), NGC2264-C which has a size of 
0.3 pc and contains two intermediate mass Class YSOs at its 
center. The here presented statistical analysis of several MDCs 
in the same molecular complex enables to establish a core mass 
function with 9 objects in the high-mass regime, and discuss the 
statistical properties of the fragmentation of MDCs. 
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Table 4. Properties of the star formation activity in the MDCs, 
and predictions assuming a normal IMF with SFE=30 %. 



Name 


Mass 
[Mo] 


^PdBI 

[Mo(%)] 


[Mo(%)] 


^MYSO 

[Mo(%)] 


^30% 

[Mo] 




[Mo] 


[Mo] 


CygX-N3 


84.0 


7.2(8.6) 


18.7(22.3) 


14.4(17.1) 


25.2 


66 


2.9 


8.1 


CygX-N12 


86.0 


17.3(20.1) 


28.8(33.5) 


24.0(27.9) 


25.8 


68 


3.0 


13.4 


CygX-N48 


197.0 


26.1(13.3) 


28.2(14.3) 


22.0(11.2) 


59.1 


156 


5.6 


10.1 


CygX-N53 


85.0 


28.5(33.6) 


50.9(59.9) 


37.2(43.8) 


25.5 


67 


2.9 


22.8 


CygX-N63 


58.0 


35.5(61.2) 


64.3(111) 


64.3(111) 


17.4 


46 


2.2 


55.2 


Total 


510.0 


114.6(22.5) 


190.9(37.4) 


161.9(31.7) 


153.0 


403 


11.3 


55.2 



^ It includes all stars and brown dwarfs from 0.01 to 120 Moamong 
which 63 % are stars. 



4. 7. Do massive fragments result from a turbulence 
regulated evolution? 

The present observations clearly demonstrate that most of the 
IMDCs are sub-fragmented and should therefore be classified 
as proto-clusters while they were found to have exactly the re- 
quired masses, surface densities, and line widths to be single pre- 
cursors of OB stars in the turbulence-regulated view of McKee 
|& Tan| pOQ2| ) (see Sect. |4.1| ). This would be acceptable if thes 
IMDCs would, at the end, contain only low-mass protostars. But 
this is not the case, these proto-clusters actually contain massive 
fragments, which have the potentiality to form OB stars. These 
fragments have the same separations, and probably roughly the 
same sizes than in nearby, low-mass proto-clusters, making them 
significantly more compact and dense than the [IMcKee & Tan 
(|2002|) cores. Such smaller sizes seem difl&cult to reconcile with 
the slow turbulence regulated evolution to turbulence supported 
cores which have to be in equilibrium with their surroundings. In 
this scenario, the size of the cores results from the Larson scaling 
laws of turbulent clouds. 

Only the most massive fragment, CygX-N63 ]V[]V[1, seems to 
be a good candidate to be a turbulence regulated massive core. 
Interestingly enough, it is also the most massive single object 
of the sample. With a size of 0.06 pc (FWHIM) and therefore a 
half-mass radius of ~ 0.03pc, it has a size very close to the core 
radius predicted in IMcKee & Tan ( 2003 ). On the other hand, this 
may just indicate that CygX-N63 is following the classical scal- 
ing laws and did not fragment further for some unknown reason, 
or that it will sub-fragment in the near future. Also some fur- 
ther investigations are required to establish in which kinematic 
state this core is. In fact, the line width in CygX-N63 is actually 
rather small for such a massive core in the quasi- static evolu- 
tion view of IMcKee & Tan| ( 2002 ). It leads to an exceptionally 
large IM/lMvir ratio of 2.7. The line width is only 2.7 times larger 
than the thermal width at 20 K. To understand the origin of such 
a core, it is therefore necessary to further investigate its proper- 
ties, and especially its kinematics, and the nature of the observed 
velocity dispersion to establish, for instance, whether or not it is 
a true isotropic turbulent dispersion. 



4.8. IMF and the core mass function of the fragments 

In the simplified view that the IMDCs are proto-clusters, they 
should fragment into a population of proto-stellar objects with 
a mass function which mimics the IIMF of stars. With a typical 
SFE of 30 %, a 100 Mq MDC would form of the order of 80 stars 
(including ~ 30 brown dwarfs) with a maximum stellar mass of 
~ 3.3 IMo for a Kroupa IMF ( |Kroupa|2QQT] ). The total number of 
stars A^^ and the mass of the most massive star Mf^^ expected to 



be formed in each IMDC, and globally for the 5 IMDCs, are given 
in Tab.|4] For comparison, the observed maximum masses of the 
fragments in each IMDC are also listed. It is clear that they are 
all systematically larger than what the IMF would predict. The 
mass function of the 23 detected fragments in the all 5 IMDCs in 
Tab.|3]is plotted in Fig.[5]in comparison with the expected global 
IIMF for the 5 IMDCs. Our PdBI data cannot detect any low-mass 
fragments, and the completeness limit is high, probably close to 
3 to 4 IM© for 4 of the 5 IMDCs, and even higher for CygX-N63 
for which it would be difficult to recognize a fragment of less 
than ~ IOMq (CygX-N63 MM2 with a mass of 9.1 Mq is ac- 
tually a highly uncertain detection). This figure illustrates that 
there is a global excess of massive fragments compared to nor- 
mal low-mass proto-clusters which mimic the IIMF. Even in the 
very un-probable case that the fragments would actually not be 
larger than the sizes derived with Gaussclump on the interfero- 
metric maps, the mass function of MpdBi is still in excess to the 
IIMF (dashed histogram in Fig.|5]). 

Globally for the 5 IMDCs, with a SFE of 30 %, only one OB 
star should statistically be formed with a mass of 11.3 Mq. The 
second most massive star should be a 6.8 Mq star. It seems there- 
fore clear that the most massive IR-quiet IMDCs of Cygnus X 
are the location of the formation of a large fraction of high-mass 
proto-stellar objects which are likely candidates to form OB stars 
in excess to a normal IlMF-like mass function. We further dis- 
cuss in the following sections the important implications of this 
result. 



4.9. Do gravo-turbulent fragmentation plus competitive 
accretion can explain the observations? 

Due to the local, low values of the Jeans mass, the IMDCs are 
actually expected to fragment and form clusters of protostars 
( [Dobbs et al. POOS ') which may later col lect more mass in th e 
central regions by competitive accretion (Bonnell & Ba te|2006| ). 
In the case of strong radiative warming from already formed 
massive stars, the local Jeans mass could be larger and the level 
of fragmentation (number of fragments) could be reduce d and 
could therefore form higher mass fragments (e.g. Krumholz &| 
Bonnell| f2QQ7 ). The local temperatures are however not high 
enough (yet?) to have Jeans masses of the order of lOlVLo (see 
Tab. [2]). Note also that the turbulent fragmentation has difl&culties 
in forming high-mass fragments leading to core mass functions 
with a deficit of high-mass cores (see discussions in Bonnell 
et al.|2QQ7 , and references therein). An additional process, such 
as competitive accretion, is required to explain high-mass stars, 
and to fully explain the observed IMF in this view. 
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The observed fragmentation of the MDCs could be seen as 
a support for the gravo-turbulent fragmentation scenario. The 
fraction of mass in the high-mass fragments, and the level of 
fragmentation (number of fragments) in the MDCs are however 
clearly not compatible with the pure gravo-turbulent fragmen- 
tation which would lead to a whole cluster of mostly low-mass 
objects (see previous section). The most striking examples are 
CygX-N12, CygX-N53 and CygX-N63 which have 28, 44 and 
100 % of their total mass, respectively, in their two most mas- 
sive fragments (see Tab. |4]) while only a few percents would be 
expected. In the context of the Bonnell and coll. view, a large 
fraction of the Cygnus X MDCs would already be in the stage 
after an intense competitive accretion to explain such high-mass 
objects. This might be surprising since the cores are still cold 
and therefore young. On the other hand, the high fragmentation 
efficiency (fraction of mass in dense fragments; see next section) 
indicates that the proto-clusters are already dominated by gravity 
of the proto-stellar objects as required for competitive accretion. 
But it would have time to work only if the crossing time in the in- 
ner region is shorter than the evolutionary time. If the evolution 
is dominated by the gravity, this evolutionary time is close to the 
global free-fall time which is typically smaller than the crossing 
time at the scale of the cores (see Sect. |4.1| ). So to have a sig- 
nificant competitive accretion, the crossing time should be sig- 
nificantly smaller in the inner regions with typically an increase 
or at least a conservation of the velocity dispersion at the small- 
est scales. This has to be further investigated. Alternatively, the 
bias of the fragmentation towards high-mass fragments could be 
primordial if, for instance, another support such as the magnetic 
field could prevent from a high degree of fragmentation as it was 
already obtained using MHD numerical codes for the collapse of 
low-mass pre-stellar cores ( Hennebelle & Teyssier|2008) . 



4.10. High fragmentation efficiencies in MDCs 

To further discuss the resulting fragmentation in the MDCs, we 
define the fragmentation efficiency as the fraction of the total 
mass in the MDCs which is in self-gravitating fragments. Since 
the targeted MDCs are cold and young, we expect that most 
of the mass of the formed stars is still in the protostellar en- 
velopes/fragments and that the fragmentation efficiency can then 
be seen as a proxy for the SFE in the MDCs. The degree or level 
of fragmentation discussed in the previous section is more re- 
lated to the number of fragments, and could be referred to the 
multiplicity of the fragmentation. It is interesting to note that 
the fractions of the single dish emissions recovered by the inter- 
ferometer are exceptionally high in our sample, and especially 
for CygX-N12, CygX-N53 and CygX-N63. This indicates that a 
large fraction of their mass is in small-scale structures. Since the 
small scale structures are expected to consist of self-gravitating 
fragments, we believe that this indicates a high fragmentation 
efficiency. To quantify this eff'ect, we have measured the frac- 
tion of the single-dish fluxes recovered by the interferometer by 
integrating the emission above 3 cr associated for all sources de- 
tected at a 5 (T level with Gaussclump in the 1.3 mm maps. The 
obtained values, expressed in terms of masses (Mf^^^) and re- 
lated fractions are listed in Tab. |4] They range between 10 and 
20 % for 3 of the 5 MDCs which are rather typical values for 
such massive cores but reach the exceptional values of 34 and 
61 % for CygX-N53 and CygX-N63. Alternatively, the fragmen- 
tation efficiency is simply evaluated as mII^^ by summing up the 
masses of the fragments listed in Tab.[3](Mtot)- This second esti- 
mate mI^^ is more indirect since it is based on the extrapolated 
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Fig. 6. Fragmentation efficiency as a function of the average den- 
sities for dense cores of similar sizes (FWHM ~ 0.1 pc) but 
of difl'erent masses in Cygnus MDCs (filled symbols) and in 
p Ophiuchus (open triangles; |Motte et aL] [l998 ). A transition 
from low efficiency cores to possible, single collapsing cores is 
observed from a few 10^ to a few 10^ cm"^ . The displayed ranges 
(uncertainty bars) for the Cygnus cores correspond to the range 
between Mf^^^ and M^^^ from Tab. Ul To help the clarity of the 
plot, the displayed location of CygX-N3 corresponds to its aver- 
age density shifted by a factor of 1 . 1 . 



masses Mtot but it may approach closer the true fragmentation ef- 
ficiencies as it is corrected for the spatial filtering of the interfer- 
ometer. We see that for CygX-N3 and CygX-N46, the efficiency 
is below 25 % and is therefore not exceptionally high. CygX- 
N12 is intermediate with an efficiency which migh be close to 
30 %. In contrast, CygX-N53 and CygX-N63 are clearly above 
30 %, indicating that the mass of these 2 MDCs is dominated 
by dense fragments, and that they may convert a large fraction 
of their mass into stars. The case of CygX-N63 can be seen, in 
this context, as the extreme case where the fragmentation was 
so efficient that all the mass is in a single centrally concentrated 
objecj^ But interestingly enough, the MDCs with the highest 
SFEs are also the ones for which most of the mass is in a few 
massive fragments with virtually 100 % of the mass in a single 
massive object for CygX-N63, 44 % of the total mass of CygX- 
N53 in only two high-mass fragments (27 % in MMl and 17 % 
in MM2), and still 28 % for CygX-N12 in two fragments. We 
have listed in Tab. |4] as M^^^^ the total mass in the high-mass 
fragments (the 9 precursors of OB stars discussed in Sect. |4.5| ) 
for each MDC. 

Also we note that the most extreme MDC, CygX-N63, has 
the highest average density (see Tab. [2]). [Motte et al^ ([T998 ) ob- 
tained, in p Ophiuchus, a similar trend of an increase of fragmen- 
tation efficiency as a function of density (see their Fig. 7). The 
fragmentation efficiencies of the MDCs are displayed in Fi g.j6| as 
a function of the average densities derived in Rgo from Tab. [2| For 
comparison, th e values obta i ned fo r similar size, but less massive 
dense cores by |Motte et al.| ( p^98] ) are plotted as open triangles. 
We have re-evaluated the densities to the way we derived them 
in the present paper for a proper comparison. It is striking that 
the Cygnus MDCs are clearly confirming the trend suggested by 



^ The obtained mI^^ = 64.3 M© is larger t han the 58 Mq obta ined 
from the gaussian fit of the MAMBO data in [Motte et aLlpOOT) ), but 
only slightly, and is therefore compatible with all the mass of the MDC 
in a single object. 
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the efficiencies in Oph A and B2 that the highest density cores 
have a larger fraction of their masses inside pre-stellar and pro- 
tostellar objects. The origin of these high efficiencies has to be 
further investigated. We can however note that the higher the av- 
erage density, the shorter the free-fall time is with the possible 
consequence that a global collapse at the scale of the core is ex- 
pected for the highest density cores. In this simplified view, the 
clear trend in Fig. [6] could correspond to the transition from low 
density cores in the standard, large scale equilibrium between 
turbulent, magnetic pressures and gravity and dense cores domi- 
nated by their self-gravity leading to a virtually 100 % efficiency 
in single, collapsing protostars. 

4.11. Primordial mass segregation in proto-clusters? 

No matter what is the origin of the most massive fragments in the 
MDCs of Cygnus X, it is clear that the population of protostars 
in the center of the observed MDCs is biased towards high-mass 
objects. This is at least certainly true for CygX-N12, CygX-N53 
and CygX-N63. It is a representative result since it was obtained 
from a complete sample of dense cores for the whole Cygnus 
X region. This clear excess of high-mass protostars in the cen- 
ter of a few MDCs in Cygnus X is actually what is expected to 
explain the primordial mass segregation of stellar clusters. The 
young stellar clusters, such as the Orion Nebula Cluster, tend to 
have their highest mass members only in their central, densest 
regions ( [Hillenbrand & Hartmann||1998| ). For the youngest and 
not too rich clusters, the observed mass segregation can not be 
due to early dynamical relaxation and is understood as the ev- 
idence for a primordial mass segregation originating from the 
star/cluster formation process itself (see also |Bonnell et al.|2007 
and references therein). Our new results on the MDCs of Cygnus 
X therefore demonstrate the existence of such a primordial mass 
segregation, at least in some MDCs/clusters. 

If the MDCs correspond to the mass segregated, inner re- 
gions of the proto-clusters, where and when do the low-mass 
stars form? Inside the most extreme MDCs of Cygnus X (CygX- 
N12, CygX-N53 and CygX-N63) which are clearly mass segre- 
gated, it seems that there is not enough mass left to form a large 
population of low-mass stars to compensate for the observed ex- 
cess of high-mass stars. So some surrounding gas has to be col- 
lected in the near future from some lower density gas around 
the MDCs to form the low-mass counterparts. We note how- 
ever that this star-forming gas cannot be located too far from the 
center of the MDCs in order that the formed stars can incorpo- 
rate the final cluster. To be dynamically linked to the cluster, we 
take as a quantitative criterium that they should have a crossing 
time several times smaller than the age of the youngest embed- 
ded clusters. With a typical velocity dispersion of 0.9 km.s"\ 
and a crossing time of 1/3 of 1 Myr, the obtained distance is 
~ 0.3 pc. The observed MDCs of Cygnus X are all embedded in 
larger scale clumps which have total masses ranging from 260 
to 7300 Mo , sizes ranging from 0.25 to 0.74 pc (FWHM), and 
average densities ranging fr om 2.5 x 10^ to 1.3 x 10^ cm"^ (see 
Tab. 2 in [Motte et al.||2007| ). It is therefore of the order of sev- 
eral hundreds of solar masses of gas at relatively high density 
which is available in a radius of ~ 0.3 pc around each MDC 
which should be enough to form a large population of low-mass 
stars in their surroundings. Alternatively, a large, already formed 
population of low-mass stars and protostars might be present in 
the immediate surroundings of the MDCs. The low-mass stars 
and protostars are difficult to observe at the distance of Cygnus 
X and it is therefore difficult to establish if there is already such 
a large enough population of low-mass counterparts in the sur- 



roundings. We will investigate this second possibility in a more 
complete way in a forthcoming paper using the recent Spitzer 
data from the legacy program towards Cygnus X ( jHora et aL 
[2009] ). 

5. Conclusions 

We have been imaging at high spatial resolution with the IRAM 
PdBI interferometer the most massive IR-quiet dense cores dis- 
covered by [Motte et al.| ( j2007| ) in Cygnus X. These MDCs are 
among the best targets to investigate the origin, and the earliest 
phases of the formation of high-mass stars and of OB clusters. 
From this unique dataset, we could reach the following conclu- 
sions and results: 

1. All except one MDC (CygX-N40) were found to be 
bright and easily detected with the PdBI in a beam of ~ 
3.0''x 2.3''at 3.5 mm, and ~ l.r'x 0.9''at 1.3 mm. 

2. There is a striking morphological diff'erence between the 
lower spatial resolution maps at 3.5 mm and the high res- 
olution maps at 1.3 mm for 3 out of the 5 MDCs. While at 
3.5 mm, the MDCs are dominated by bright central emis- 
sions, at 1.3 mm, these emissions clearly split into a pop- 
ulation of compact fragments with typical separations of 
5000 AU and sizes larger than 1000 AU. 

3. All except one MDC (CygX-N63) are sub-fragmented. 
CygX-N63 consists of a single, centrally concentrated ob- 
ject down to the smallest scales. It could actually be a single 
massive protostar with an envelope mass as large as ~ 60 M© 
and may therefore form an O star. 

4. A population of 23 compact dense fragments are recognized 
as certainly self-gravitating, proto-stellar objects. Among 
them, 9 objects are good candidates to be true, individual 
high-mass protostars in the sense that they are precursors of 
OB stars. 

5. The targeted MDCs of Cygnus X have gl obally the prop- 
erties proposed as initial conditions by McKee & Tan 
( 2003 ) for the turbulence-regulated, monolithic collapse. But 
in contrast to this scenario, these MDCs are mostly sub- 
fragmented. Only CygX-N63 is found to be single. To further 
investigate whether this particular core could have evolved 
towards global collapse by a slow evolution, regulated by a 
true turbulent field, a detailed study at high spatial resolution 
is mandatory to, for instance, probe its kinematical structure, 
and the nature of its velocity dispersion. 

6. On the other hand, the properties of the fragmentation of the 
MDCs are not in agreement with a pure gravo-turbulent frag- 
mentation. A large fraction of the mass of 3 out 5 MDCs is 
located in a few, high-mass objects which would not be ex- 
pected in the Jeans mass regulated fragmentation. The local 
Jeans masses are too low, of the order of 1 M© . Similarly, 
the core mass function of the fragments is globally (for the 
5 MDCs) biased towards high-masses as compared with a 
normal IMF-like core mass function. 

7. We found a clear trend that the fragmentation efficiency, de- 
fined as the total fraction of mass in self-gravitating frag- 
ments, is increasing with the average density of the cores. 
Since the density is directly related to self-gravity, we inter- 
pret this result as due to the larger eff'ect of self-gravity in the 
densest cores leading to global collapse of the most extreme 
cases. 

8. In the same time, as the density increases, the fragmentation 
tends to lead to a fewer number of predominantly high-mass 
fragments in the central part of the cores. The central parts 



14 



S. Bontemps, F. Motte, T. Csengeri, N. Schneider: Fragmentation in the massive dense cores of Cygnus X 



of at least 3 ot the 5 MDCs are so clearly dominated by mas- 
sive objects that they are certainly leading to regions with a 
clear mass segregation. This mass segregation is proposed to 
correspond to the expected primordial mass segregation of 
young stellar clusters. 
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